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Quantum state measurements using multi-pixel photon detectors 
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The characterization and conditional preparation of multi-photon quantum states requires the use 
of photon number resolving detectors. We study the use of detectors based on multiple avalanche 
photodiode pixels in this context. We develop a general model that provides the positive operator 
value measures for these detectors. The model incorporates the effect of cross-talk between pixels 
which is unique to these devices. We validate the model by measuring coherent state photon number 
distributions and reconstructing them with high precision. Finally, we evaluate the suitability of 
such detectors for quantum state tomography and entanglement-based quantum state preparation, 
highlighting the effects of dark counts and cross-talk between pixels. 

PACS numbers: 42.50.Dv,42.50.Ar 



I. INTRODUCTION 

The development of detectors with photon number 
resolution in the 'few' photon regime is of key interest 
in quantum optics. Such detectors play a central role 
in the preparation and characterization of non-classical 
multi-photon states. Experiments requiring photon num- 
ber resolution include the generation of entangled pho- 
ton sources via spontaneous parametric down conver- 
sion (SPDC) with strong pump fields [![, linear optics 
quantum computation [2| and a proposal for beating the 
Heisenberg limit of intcrferometry with mixed SPDC and 
coherent state fields 0. 

The realization of photon number resolving detectors 
is extremely challenging. The benchmark properties for 
these are quantum efficiency, dark count rate, dead time, 
operating temperature and the degree and quality of pho- 
ton number resolution [J] . One of the first realizations of 
a photon number resolving detector was the 'visible light 
photon counter' (VLPC) Q, an avalanche photo-diode 
whose output pulse height is proportional to the num- 
ber of detected photons due to low multiplication noise 
in the avalanche process. It has extremely high detec- 
tion efficiency of around 90% but requires cryogenic cool- 
ing and exhibits a high dark rate. More recently, super- 
conducting devices have been investigated. These devices 
have high quantum efficiencies and low dark-counts, but 
they too require cryogenic coolingto reach extremely low 
working temperatures [f| 0, B In addition, a few 
groups have realized a limited degree of photon number 
resolution using standard detectors at room temperature 
[ToL [TH . Another approach involves the use of indepen- 
dent single photon detectors at the output ports of a 
balanced array of beam-splitters [HI, [l3| . Although the- 
oretically tractable, this solution is experimentally unap- 
pealing due to the large number of detectors involved and 
the complexity of the setup. A practical way to overcome 
this difficulty has been the use of temporal modes instead 
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of spatial ones [lj, [l5| ■ One can separate photons into 
16 separate temporal modes using only two detectors. 

An alternative method for separating photons into dis- 
tinct spatial modes is allowing the beam to spread via free 
space propagation onto an array of single photon detec- 
tors. Each detector in the array provides a binary 'click' 
when one or more photons are detected by it. When the 
total number of photons is significantly smaller than the 
number of detectors then it may be assumed that each 
pixel detected at most one photon. In this regime, the 
number of pixels that 'clicked' is equal to the number 
of detected photons. Clearly, high efficiency avalanche 
photo-diodes (APDs) are favorable candidates for pixels. 
Recently, the fabrication of an array of APD pixels for 
multi-photon detection has been realized [H, [TtJ ■ These 
detectors which we refer to as 'multi-pixel detectors' have 
been studied extensively for use in particle physics experi- 
ments such as the International Linear Collider (ILC) and 
T2K [l8j . In this paper we examine, both experimen- 
tally and theoretically, the use of multi-pixel detectors 
from the point of view of the quantum optics commu- 
nity. These detectors are appealing in that they exhibit 
a high single pixel efficiency, fast detection rates and ex- 
cellent photon number resolution. Two limitations which 
must be taken into consideration are the relatively high 
dark-count rate and the possibility of cross-talk between 
pixels. We note that an intensified CCD camera has been 
previously used for photon counting with number reso- 
lution 191. Such cameras however, have slow repetition 
rates and a quantum efficiency lower than APDs. 

In sections |TT] and IIIII we describe the experimental 
setup and our method for data acquisition which min- 
imizes the room temperature dark-count rate. In section 
IIVI we develop a theoretical model for relating the mea- 
sured data to the actual photon number distribution. In 
section [V] we validate the model by measuring coherent 
state statistics. Finally, in section IVI1 we evaluate the 
applicability of multi-pixel detectors to quantum tomog- 
raphy and entanglement-based quantum state prepara- 
tion. 
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II. EXPERIMENTAL SETUP 

The Hamamatsu multi-pixel photon counter (MPPC) 
serves as our prototype detector. We use a module 
(C10507-11-050U) that incorporates the 400 pixel MPPC 
and peripheral electronics. Our experimental setup is 
shown in Fig. [TJ The light source is a mode locked 
Ti-Sapphire oscillator (Spectra Physics, Tsunami) which 
emits 120 fs pulses centered at 810 nm with a 80 MHz 
repetition rate. The detector's output pulses are ~ 20 ns 
FWHM, therefore the maximum usable laser rep. rate 
is ~ 10 MHz. We use a pulse picker (PulseSelect, APE- 
Berlin) to reduce the laser repetition rate to 1 MHz. Af- 
ter strongly attenuating the beam using ND filters we al- 
low the remaining light to impinge on a multi-mode fiber 

wtiirln iq rnnnlprl HirApt.lv tn tVip r\c*tprtnr TtiA A&t&rtrtr 
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FIG. 1: Setup for evaluation of MPPC. 



FIG. 2: Examples of waveforms obtained from the MPPC. 
The vertical line marks the arrival time of the illuminating 
pulse, (i) Waveforms which are rejected since they are overrun 
by dark-counts, (ii) Selected waveforms used in the data- 
analysis: (a) A two avalanche waveform; (b) A one avalanche 
waveform; (c) A zero avalanche waveform which starts to rise 
outside the designated time-window due to a dark-count; (d) 
A zero avalanche waveform that is flat during the whole time 
window; 



III. DATA AQUISITION 

The output of a multi-pixel detector is the sum of 
the outputs of each of the pixels. Thus the height of a 
given output pulse is proportional to the number of pix- 
els which 'clicked' in a given detection event. Multi-pixel 
detectors typically suffer from a high overall dark count 
rate. This is because the total dark count rate is a sum of 
the dark counts of all the pixels. In our experiment the 
computer based digitizer receives a gating signal from the 
pulse picker which is synchronized with the pulse arrival 
times. For each trigger we acquire a 20 ns waveform at a 
1GHz sampling rate in pre-trigger mode starting from 5 
ns before the trigger arrival time. We use the following 
scheme to post-select valid waveforms: 

1. If the waveform is not at the digitizers' zero level 1 
ns before the arrival of the trigger then the signal is 
discarded. This implies that a dark-count occurred 
prior to the trigger and may interfere with the read- 
out process by obscuring the signal from the actual 
pulse or causing an undesirable after-pulse (see Fig. 
[Hi)). Otherwise continue to step two. 

2. If a rising edge occurs within 3 ns of the trigger 
arrival then the waveform value is saved at some 
time close to its peak (e.g. t = 5 ns). This value is 
later used for determining the number of photons 
in the pulse. If no rising edge occurs then the value 
is saved, corresponding to a zero photon event 
since no pixels fired in the designated time window 
(see for example Fig. [U(ii) curves c and d). 



This detection scheme ensures the rejection of most 
of the dark-count and after-pulse events. The obtained 
single photon dark count rate is 2.3 x 10~ 3 counts per 
pulse. The dark count rate without the post-selection 
scheme described above is 2.2 x 10~ 2 i.e. almost an order 
of magnitude higher. We note that dark counts ma y b e 
significantly reduced using liquid nitrogen cooling [2fj |. 
Using the described method we process data at a rate of 
80 KHz, limited mainly by the digitizers' internal buffer 
size. Using the largest available buffer available for the 
NI PCI-5152 digitizer series would allow a much higher 
data processing rate. 

Fig. [3] shows a histogram of the pulse heights obtained 
when illuminating with a coherent state. The MPPC is 
shown to exhibit a high degree of photon number reso- 
lution manifested in the well separated peaks. Events of 
up to 10 pixels firing simultaneously are easily resolved. 



IV. THEORY AND DATA-ANALYSIS 

The probability p' n for a given number of avalanches 
per pulse is proportional to the area underneath the re- 
spective peak in the histogram in Fig. [3j In this sec- 
tion we solve the problem of finding the relation between 
p' n and the actual photon number distribution imping- 
ing on the detector which we denote p n . This relation is 
a linear transformation and can be conveniently written 
in matrix form [3, EH- F° r multi-pixel detectors, the 
transformation has three contributions originating from 
losses, cross-talk and dark-counts. 
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FIG. 3: A pulse height histogram obtained using a coher- 
ent state light source. The well distinguished peaks indicate 
accurate photon number resolution. 



by s D , the relation for dark-count compensation is 

p' = M D p 

!1 — e D if m = n; 
e D if m = n - 1; 

otherwise. 

n,m = 0,1, 2... (2) 

Note that e D is the dark-count probability which would 
occur in the absence of cross-talk (see Eq. J?])). 



C. Cross-talk 

Due to the proximity of pixels, an avalanche in a given 
pixel may induce avalanches in its neighbors. We use a 
model with one parameter, e XT , which is the probability 
that a given avalanche will induce cross-talk. We find the 
following relation for cross-talk compensation: 



A. Loss 

The photon detection efficiency (PDE) of a multi-pixel 
detector can be obtained by multiplying the PDE of a 
single pixel by a geometric factor that accounts for the 
unutilizeduneffcctive area between the pixels. According 
to specifications [ll[, the geometric factor for the 400 
pixel device is 61.5% and the total PDE, r\, has a peak 
value of 50% at 400 nm and drops to 8% at 800 nm. 
Using the definitions introduced above, we find that loss 
can be compensated using the following vectorial relation 
between p and p': 

p' = M L p 

n,m = 0,1,2..., (1) 

where Ml is a square matrix used to describe the effects 
of loss. 



p' = M X T P 

[MxT] n , m =( n ! m ) e n x - m (l-e XT ) m 

n,m = 0,1,2... . (3) 

Here, Mxt is the cross-talk matrix. Note the similar- 
ity of the matrix elements of Mxt and Ml . The main 
difference being that Mxt is upper trigonal and Ml 
is lower trigonal implying that loss transfers probability 
from high photon numbers to lower ones whereas cross- 
talk does the opposite. The com pen sation for cross-talk 
was dealt with previously by fl6l.l20l| in the form of a re- 
cursion relation. Our formalism allows for more than one 
cross-talk event per-pulse. For example, three photons 
may initiate up to six avalanches. The matrix for first- 
order-only cross-talk may be constructed by using the 
diagonal and sub-diagonal elements of Mxt only. We 
show below (Fig. that high order cross-talk events are 
crucial for fitting high photon numbers whose probabili- 
ties are low. In addition, the matrix form of our method 
allows the relation between p and p' to be conveniently 
inverted (see Section HVD[) . 



D. Photon Number Reconstruction 



B. Dark Count 

Dark counts are avalanches in APD pixels that origi- 
nate from thermal noise. In principle a dark count can 
consist of multiple avalanches from different pixels. We 
consider only single avalanche dark counts since the effect 
of higher order dark counts is negligible. Denoting the 
probability per pulse for a single avalanche dark count 



There are two approaches to reconstruction of the pho- 
ton number distribution. The first approach requires as- 
sumed knowledge about the photon number distribution 
of the light illuminating the detector (i.e. coherent state, 
spontaneous parametric down-conversion etc.). Usually 
the distribution will have a parameter, such as mean pho- 
ton number, which must be chosen correctly in order to 
fit the measured data. The relation between the assumed 
distribution p and the measured data p' is given by ap- 
plication of the loss, dark-count and cross-talk matrices 
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in the correct order, 



p' = M XT M D M L p. 



(4) 



Using Eq.(UJ) we verify that the assumed distribution is 
compatible with the measured data thus allowing us to 
find the distribution's parameters. 

The second approach for retrieving the photon number 
distribution requires no a-priori assumptions. In this case 
we use the inverse of the relation in Eq. (0| to obtain p 
from p', 



p = M^ 1 M D ' ■ M x l T ■ p' . 



(5) 



All of the matrices involved can be inverted analytically 
as shown previously for the case of loss and dark counts 
[22| . This second approach is less stable than the first, 
since the inverted relation in Eq. ([5]) is very sensitive 
to noise in the measured statistics. For photon numbers 
with small probabilities this can yield un-physical oscil- 
lations in the obtained photon number distributions and 
negative probabilities. The inversion may be stabilized 
using physical assumptions about the result [23|, [24[ • 



number probabilities of the incoming state by p. Us- 
ing Eq. (|3]) we find the following relation between the 
first two elements of these vectors, 



Vx = (pi (1 - Er>) + PO £n J ( 

Po = Po (1- £d) ■ 



I -By 



(8a) 
(8b) 



We assume that p = exp(— (n)) and pi = cxp(— (n)) x 
(n), as expected from a coherent state. As a result, 
we find that Eqs. (I8a|8b|) form a set of two equations 
with two unknowns, (n) and e XT and can be solved to 
determine both of them. Using this method we cal- 
culated e XT for a number of different values of (n) in 
the range 0.5 — 3.5. As shown in Table [V] we obtained 
Ext = 0.0975 ± 0.0015. This method for obtaining e XT 
is more accurate than using only two photon dark-count 
rate as a measure of cross-talk 171. 



(n) 0.86 


1.66 


2.30 


3.13 


e XT 0.0960 


0.0985 


0.0965 


0.0965 



TABLE I: Measured cross-talk probability, obtained for co- 
herent states with different mean photon numbers. 



V. COHERENT STATE STATISTICS AND 
DETERMINING THE CROSS-TALK 
PROBABILITY 

It is well known that the state of light emitted by a 
pulsed laser is a coherent state with the following photon- 
number statistics, 



Pk = e-<"> 



(n) k 
k\ 



(6) 



where (n) is the average photon number. The applica- 
tion of loss transforms one coherent state to another with 
the mean photon number multiplied by the efficiency r). 
It is therefore impossible to calibrate the efficiency of a 
detector using coherent states without the use of a ref- 
erence detector [23|. We do not, therefore, account for 
the effects of loss in this measurement. Cross-talk, on the 
other hand, changes the form of a coherent state distribu- 
tion and allows us to accurately determine the cross-talk 
probability by measuring coherent states. To do this, we 
first take a measurement of the dark counts by block- 
ing the laser in our setup. Denoting the measured single 
avalanche dark-count probability by e' D it can be shown 
that the dark count rate deducting the effect of cross-talk 
is 



(7) 
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Next, we take a measurement of the laser pulses. We 
measured 10 6 pulses in ~ 10 sec. Following the nota- 
tion of the previous section, we denote the measured 
avalanche number probabilities by p' and the photon 
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FIG. 4: Measurement of a coherent state, (a) The raw 
avalanche statistics, (b) An uncorrected coherent state with 
(n) = 1.66 (c) The data in (b) multiplied by the matrix 
M X T ■ M D with first order cross-talk compensation only, (d) 
The same as (c) but with full cross-talk compensation. 



The data measured for a coherent state with (n) = 1.66 
together with the compensated theoretical distribution is 
shown in Fig. [4] There is very close agreement between 
the experimental data and the theory over more than four 
and a half decades of probability values. The theoretical 
distribution with compensation for only first order cross- 
talk (i.e. up to one cross-talk event per pulse) is shown 
to illustrate that high order cross-talk compensation is 
necessary for obtaining an accurate fit at large photon 
numbers. 
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VI. QUANTUM STATES: TOMOGRAPHY AND 
CONDITIONAL PREPARATION 

The complete characterization of quantum states of is 
known as quantum state tomography. Typically, when 
characterizing a state of light, it is assumed that there is 
an unlimited amount of copies at our disposal. A large 
number of measurement results are used to determine 
the most likely input state. In the following we use 
the term click when referring to the discrete output of 
a general photon number resolving detector. An ideal 
detector would give exactly n clicks when illuminated 
by an n photon Fock-state. Unfortunately, existing de- 
tectors have various properties which obscure the one- 
to-one correspondence between the number of clicks and 
the number of impinging photons. The main properties 
that determine a given detector's applicability for quan- 
tum tomography are single photon detection efficiency 
and dark count rate. For time or space multiplexed de- 
tectors, the number of different detection modes is very 
significant. Since each mode can detect at most one pho- 
ton, a large number of modes ensures that two photons 
will almost never reach the same mode. The large num- 
ber of pixels in the MPPC's array make it very appeal- 
ing in this respect. The existence of cross-talk, which is 
unique to multi-pixel detectors, must also be accounted 
for. The complete characterization of a quantum detector 
isgiven by its positive operator value measures (POVM) 
|23j |. Given an input density matrix p, the probability 
Pn iP of obtaining a detection outcome n is 



Pn, P = tr [p7r n ] 



(9) 



where {7r„} is the POVM. For a phase independent de- 
tector, 7r„ is diagonal in the Fock state basis and has the 
general form, 



fc=0 



( k n) \k){k\ 



(10) 



Here, 0^ is the probability to obtain n clicks when the 
detector is illuminated by a k photon Fock state. It was 
shown in Eq.(T4]) that the transformation between the il- 
luminating state and the click statistics is given by, 




12 16 
Fock State I 



FIG. 5: Elements of the MPPC's POVM, based on the model 
developed in Section IIVI The different curves correspond to 
events with different numbers of clicks. Results where calcu- 
lated assuming a 50% efficiency. 



to lower Fock states since an n click event may origi- 
nate from a Fock state with less than n photons. Due to 
the large number of spatial modes, the POVM's shape 
is determined mainly by the detection efficiency. For 
time-multiplexed detectors the limited number of modes 
manifests itself in a broadening of the POVM curves and 
a significant shift towards higher photon numbers [23[. 
We conclude that multi-pixel detectors are applicable to 
quantum state tomography. The existence of cross-talk 
doesn't pose a significant limitation and the large number 
of pixels enhances the ability to reconstruct high photon 
number states. 

In entanglement-based quantum state preparation a 
measurement on one mode of an entangled two mode 
state predicts the existence of a desired state in the other 
mode [12l l22j . Quantum state preparation is most com- 
monly realized using two-mode spontaneous parametric 
down-conversion (SPDC). The wave-function of the two 
mode state before detection is given by '25] 



M 



TOT 



= M X t • M D ■ M 



L ■ 



(11) 



The matrix Mtot contains all the information needed 
for constructing the MPPC's POVM. The connection be- 
tween the two being given by, 



\M> 



TOT \k ■ 



Jn) 
7 k ■ 



(12) 



The elements of the POVM with an assumed efficiency 
•q are shown in Fig. [5] The normalization is such that 
9^ — 1. Note that in the absence of cross-talk, 



<L~tn k 



the peak of the n click curve occurs at i ~ njr\. Due 
to cross-talk, the peak is shifted slightly to the left i.e. 



W) = X! C n\n) a \n)i 



n=0 



C n = — (-l)V ne tanh n r 

coshr 

(n} a = (ri)b = sinh 2 r. 



(13) 



The detection of n photons in mode a with an ideal de- 
tector indicates the existence of n photons in mode b. 
These are so-called heralded n photon Fock states. For 
a realistic detector one defines the fidelity Q(k\k) as the 
probability that given k clicks in mode a there are actu- 
ally k photons in mode b. In terms of the POVM elements 
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k = 1 it can be seen that approaching small values of (n) 




the fidelity starts to curve downwards and finally reaches 
zero at values when (n) is comparable to the dark count 
rate. This reduction of fidelity occurs due to dark-counts 
and has nothing to do with cross-talk. Nevertheless, for 
mean photon numbers larger than O.f the MPPC obtains 
better fidelity than a single APD despite the dark-counts. 
For k = 2 the situation is even less favorable. The reason 
being that single photon events are mistaken for two pho- 
ton events due to cross-talk. We conclude that cross-talk 
significantly impedes the use of multi-pixel detectors for 
quantum state preparation at least in realistic scenario 
described above. The elimination of cross-talk as well as 
the reduction of dark counts would be desirable for this 
application. 



VII. CONCLUSION 



FIG. 6: Fidelity of the prepared state for a single photon Fock 
state (left pane) and a two photon Fock state (right pane), 
as a function of the average number of photons in either of 
the modes, calculated results. The post-selection scheme is 
described in the text. The fidelity is plotted for three dif- 
ferent heralding detectors: (a) the complete MPPC model, 
(b) the MPPC model with loss but without cross-talk and 
dark-counts , (c) a single APD (left pane) two APDs' and a 
beamsplitter (right pane). 



Q(k\k) is written as 



Q(k\k) 



Jk) 

7 k 



|C fc | 



To illustrate the effect of dark-count and cross-talk on 
the preparation fidelity, Q(k\k) is plotted in Fig [5] for 
k = 1,2. For each value of k we plot the fidelity for 
both the complete MPPC model and the MPPC without 
cross-talk and dark-counts. This can be done by omitting 
the dark count and the cross-talk matrices from M^ot- 
In the absence of a photon number resolving detector 
one could detect a single photon using a single APD and 
two photons using two APDs' and a 50/50 beamsplitter 
[2r| . The fidelity of state preparation with heralding per- 
formed using these basic solutions (assuming zero dark 
count and 50% efficiency) is plotted for comparison. For 



We have studied the use of multi-pixel detectors in the 
context of quantum optics using the Hamamatsu MPPC 
as a prototype. This device is shown to have excellent 
photon number resolution. Using triggered data acquisi- 
tion and post selection we achieved a room temperature 
dark-count rate of ~ 2 x 10~ 3 per pulse which may be 
further reduced by cooling. In addition, we experimen- 
tally measured the cross-talk probability with high preci- 
sion using coherent states. A complete model for photon 
number reconstruction was developed which includes the 
effect of high order cross-talk. Due to the large number of 
pixels, such detectors are particularly suitable for mea- 
surement of states with high photon numbers, making 
them attractive candidates for quantum state tomogra- 
phy. On the other hand, the fidelity of quantum states 
prepared using a standard entanglement-based protocol 
is found to be low due to the effects of cross-talk and 
dark counts. The full potential of these detectors will be 
realizable when cross-talk between pixels is eliminated. 
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